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Abstract. Cen A has been recently detected in the high-energy (HE) and very high energy 
(VHE) y-ray domain by Fermi-LAT and H.E.S.S. We review the observed characteristics 
and suggest a scenario where the VHE emission originates from the vicinity of the black 
hole. Motivated by the possible association of some ultra-high energy (UHE) cosmic ray 
(CR) events with Cen A, we further analyze the acceleration efficiency for a number of a 
mechanisms (including shock and stochastic acceleration), showing that most of them (apart 
perhaps from shear) have serious difficulties in accelerating protons beyond a few 10 19 eV. 
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1. Introduction 

Cen A is the nearest active galaxy, less than 
4 Mpc away. It has a peculiar radio morphol- 
ogy with two jets emerging from its nucleus, 
and giant radio lobes stretching out to 250 
kpc and extending over a 8 x 4 degree field 
on the sky (for review, see e.g. Israel 1998). 
VLBI studies have shown that Cen A is a 
non-blazar source with its jet inclined at view- 
ing angles i ~ 50° and characterized by moder- 
ate bulk flow speeds of Uj ^(0.1- 0.5) c only 
(Tingay et al. 1998; Hardcastle et al. 2003; 
Miiller et al. 201 1). Cen A harbors a supermas- 
sive black hole (BH) of mass wbh = (0.5 - 3) x 
1O 8 M (Marconi et al. 2006; Neumayer et al. 
2007). Given its estimated bolometric luminos- 
ity ~ 10 43 erg/s (Why song & Antonucci 2004), 
Cen A is rather under-luminous and believed 
to accrete at sub-Eddington rates. If its inner 
disk would remain cooling-dominated (stan- 
dard disk), accretion rates m ~ 10~ 3 /«Edd and 
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equipartition magnetic field strengths close to 
the BH of B = (2L b /r 2 g c) l/2 a 2 x 10 3 G 
are expected (r g = GM/c 2 is the gravitational 
radius). If the disk would switch to a radia- 
tively inefficient (RIAF) mode, characteristic 
field strengths may reach Bo ~ 10 4 G. 

2. Gamma-rays from Cen A 

Cen A was the only non-blazar AGN detected 
at MeV to GeV energies by all instruments on 
board the Compton Gamma-Ray Observatory 
(for review, see Steinle 2010). Fermi-LAT has 
recently reported HE (> 0.2 GeV) y-rays from 
both the giant radio lobes and the "core" (i.e., 
within ~ 0.1°) of Cen A (Abdo et al. 2010a,b): 
Both lobes have been detected up to 3 GeV, 
with step spectral slopes (photon indices close 
to 2.6), and contribute more than one-half to 
the total HE source emission. The HE lobe 
emission can be modeled as due to Compton 
up-scattering of CMB (e = 8 x 10~ 4 eV) and 
infrared extragalactic background photons by 
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electrons with Lorentz factors y e = 6 x 10 5 , 
assuming fields strengths B 0.9/l/G. This 
would imply a total energy (assuming a neg- 
ligible proton contribution) in both lobes of 
E, - 10 58 erg, and require a jet kinetic power 
Lj - 8 x 10 42 erg/s close to the available accre- 
tion power. Fermi-LAT has also reported HE 
emission up to 10 GeV from the core, again 
with a steep photon index of ~ 2.7 and with ap- 
parent (isotropic) luminosity L(> 0.1 GeV) 
4 x 10 40 erg/s. The HE light curve (using 15 d 
bins) is consistent with no variability. 
VHE (> 0.1 TeV) emission up to 5 TeV has 
been detected by H.E.S.S. in more than lOOhof 
data taken in between 2004-2008 (Aharonian 
et al. 2009). The VHE spectrum is consistent 
with a power-law of photon index 2.7 + 0.5, 
and the apparent (isotropic) luminosity is L(> 
250 GeV) ^ 2 x 10 39 erg/s. No significant vari- 
ability has been found. 

The nuclear SED of Cen A, based on non- 
simultaneous data, shows two peaks, one at 
several 10 13 Hz and one at around 0.1 MeV 
(Chiaberge et al. 2001; Meisenheimer et al. 
2007; Abdo et al. 2010b). The SED below a 
few GeV is satisfactorily described by a one- 
zone synchrotron self-Compton (SSC) model 
(e.g., Chiaberge et al. 2001). As it turns out, 
however, the same approach fails to account 
for the TeV emission observed by H.E.S.S 
(Abdo et al. 2010b). In fact, a simple extrapo- 
lation of the Fermi (power law) spectrum tends 
to under-predict the observed TeV flux. This 
could indicate an additional contribution to the 
VHE domain beyond the conventional SSC 
jet emission, emerging at the highest energies. 
Non-thermal processes in the black-hole mag- 
netosphere could offer a plausible explanation 
for this (Rieger & Aharonian 2009): Provided 
the inner disk in Cen A is radiatively ineffi- 
cient (ADAF-type), electrons can be centrifu- 
gally accelerated along rotation magnetic field 
lines to y e oc 1/(1 - r/i\) ~ 5 x 10 7 while 
approaching the light cylinder r L = c/Q, and 
thereby enable Compton up-scattering of sub- 
mm ADAF disk photons to the TeV domain, 
satisfying the observed VHE spectral and lu- 
minosity constraints. If the inner disk is of the 
ADAF-type, these VHE photons can also es- 
cape yy-absorption. Observationally, the nu- 



clear SED of Cen A peaks in the mid-infrared, 
with an apparent (isotropic) spectral luminos- 
ity of ~ 6 x 10 41 erg/s at hv ~ 0.15 eV and ev- 
idence for an exponential cut-off (!) towards 
higher frequencies (Whysong & Antonucci 
2004; Meisenheimer et al. 2007). This emis- 
sion is usually believed to be produced on 
larger scales, either by a non-thermal (non- 
isotropic!) synchrotron jet component at a dis- 
tance ~ 0.03 pc (e.g., Meisenheimer et al. 2007) 
or a (quasi-isotropic) dusty torus on scales ~ 
0.1 pc or larger (e.g., Radomski et al. 2008), 
thereby enabling sufficient dilution such that 
VHE photons are able to escape. 



3. UHE cosmic rays from Cen A 

The apparent clustering of UHECRs along 
Cen A has renewed the interest into nearby 
AGN as potential UHECR accelerators. In 
2007, the Pierre Auger (PAO) Collaboration 
initially reported evidence for an anisotropy at 
the 99% confidence level in the arrival direc- 
tions of cosmic-rays with energies ~ 6 x 10 19 
eV (Abraham et al. 2007). The anisotropy was 
measured by the fraction of arrival directions 
that were less than ~ 3° from the positions 
of nearby AGN (within 75 Mpc) from the 
VCV catalog. While this correlation has be- 
come weaker given the now available (twice as 
large) data set, the updated analysis still sug- 
gests that a region of the sky around the po- 
sition of Cen A has the largest excess of ar- 
rival directions relative to isotropic expecta- 
tions (Abreu et al. 2010). This obviously mo- 
tivates a theoretical investigation of possible 
UHECR acceleration sites in Cen A. Below 
we analyze the efficiency constraints expected 
for a number of acceleration mechanisms when 
applied to Cen A. As it turns out, most mech- 
anisms have serious difficulties in accelerating 
protons beyond a few 10 19 eV. While the ex- 
perimental situation is not fully conclusive yet, 
this result may fit into recent PAO indications 
for an increase of the average mass composi- 
tion with rising energies up to £ =* 10 19 6 eV 
(Abraham et al. 2010). 



EM. Rieger: Cen A as y- and UHE cosmic-ray Source 



3 



3.1. CR acceleration in the BH vicinity 

Rotating magnetic fields, either driven by the 
disk or the BH itself, could facilitate accelera- 
tion of charged particles: 
(i) Direct electric field acceleration: If the 

BH is embedded in a poloidal field of strength 
B p and rotating with angular frequency Q H , it 
induces an electric field of magnitude \E\ ~ 
(Q//0/)B ; ,/c. This corresponds to a voltage 
drop across the horizon r H of magnitude O ~ 
r H \E\. For Cen A, this voltage drop becomes 



<D ~ 3 x 10 19 a 



'«BH 



\1O 8 M /\1O 4 G 



[V], (1) 



where < a < 1 denotes the dimensionless 
spin parameter. If a charged particle (charge 
number Z) could fully tap this potential, accel- 
eration to E = ZeO ~ 3 x 10 19 Z eV may 
become possible. This would favor a rather 
heavy composition (e.g., irons instead of pro- 
tons) above £ c = 5x 10 19 eV. Yet, whether 
such energies can indeed be achieved seems 
questionable: First, the charge density pro- 
duced by annihilation of MeV photons emit- 
ted by an ADAF in Cen A most likely exceeds 
the Goldreich-Julian (GJ) density required to 
screen the electric field (Levinson & Rieger 
2011). A non-negligible part of the electric 
field would then be no longer available for par- 
ticle acceleration. Secondly, even if screening 
could be avoided, curvature losses (Levinson 
2000) would constrain achievable energies for 
protons to 



p \io 8 mJ 



1/2 



B„ 



\10 4 G 



3/4 



eV. (2) 



Thirdly, large-scale poloidal fields with 
strengths B p ~ 10 4 G would be required. This 
seems overly optimistic, at least for a standard 
disk (Livio et al. 1999). Fourthly, one would 
need a 1 although rather moderate spins 
are expected for FR I sources (Daly 2011). 
Therefore, efficient DC acceleration of protons 
to E c and beyond in Cen A is unlikely, but 
could be possible for heavier elements, 
(ii) Centrifugal particle acceleration: Even 
if the charge density would exceed the GJ 
density, centrifugal particle acceleration 



along rotating magnetic field lines could still 
occur (e.g., Osmanov et al. 2007; Rieger 
& Aharonian 2009). Yet, requiring that the 
acceleration timescale remains larger than the 
inverse of the relativistic gyro-frequency, CR 
Lorentz factors are limited to 

where r L is the light cylinder radius (typically 
of a few r g ). This suggests that centrifugal ac- 
celeration is unable to produce UHECRs. 



3.2. Fermi-type CR acceleration in the 
jets and beyond 

Suppose instead that CR acceleration is Fermi- 
type, i.e., due to multiple scattering off moving 
magnetic inhomogeneities, with a small energy 
change in each scattering event. We may then 
distinguish the following scenarios: 
(i) Diffusive shock (1st order Fermi): In this 
case, energetic charged particles are assumed 
to pass unaffected through a shock front and, 
by being elastically scattered in the fluid on ei- 
ther side, to cross and re-cross it several times. 
Sampling the difference Au in flow velocities 
across a shock (always head-on), the charac- 
teristic energy gain for a particle crossing the 
shock, becomes 1st order, i.e., Ae/ei cc (Au/c). 
As this is acquired during a shock crossing 
time t c ~ A/u s (with u s the shock speed and A 
the scattering mean free path), the characteris- 
tic acceleration timescale (for a non-relativistic 
shock) becomes 



^acc — 



(de/dt) 



(4) 



We can equate f acc with the timescale for 
cross-field diffusion out of the system, t e ~ 
r^,/(Ac), or the dynamical timescale, t& ~ z/u s 
(whichever is smaller), to obtain an estimate 
for the maximum achievable particle energy, 
£„ - ZeBr w B s (cf. Hillas 1984), i.e. 



^max — 2 X 10 Z 



\10 4 G/\0.1/ 



eV, 



(5) 



assuming A ~ r gyro , with r gyro the gyro-radius, 
B s = ujc, and B(z) - 4 Bq (r g /zaj) for the 
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typical magnetic field strength at location z (al- 
lowing for magnetic field compression by a 
factor of 4). Here, a, is the jet opening an- 
gle. Radio observations of Cen A indicate bulk 
flow speeds, both (!) on sub-pc and hundreds of 
pc scales, that are only mildly relativistic. This 
suggests only moderate internal shock speeds, 
B s ~ 0.2. Modest shock speeds are also sup- 
ported by the nuclear SED of Cen A with a 
synchrotron peak below 10 20 Hz (cf. Lenain 
et al. 2008) as synchrotron-limited electron 
shock acceleration results in a (magnetic field- 
independent) peak at ~ 3 x 10 19 (J3 S /0.1) 2 Hz. 
Thus, efficient shock acceleration of protons 
to energies E c and beyond seems unlikely. 
Moreover, it can be shown that otherwise also 
a jet power well in excess of the one expected 
for Cen A as an FR I-type source would be re- 
quired (Rieger 2009). 

(ii) Stochastic 2nd order Fermi: In general, 
the average energy gain due to scattering off 
randomly moving magnetic inhomogeneities 
(waves) is only second order, i.e., Ae/ei oc 
(u/c) 2 . As the energy gain is acquired over a 
mean scattering time t s ~ A/c, the characteris- 
tic acceleration timescale is 

facc ~ JdTW) ~ \v) ~c' (6) 

assuming that scattering is due to Alfven waves 
moving with u-va - B/ ^J4np. Neglecting 
radiative losses, particle energies are limited by 
escape via cross-field diffusion to 

£ <2xlO»z(-?-)pq(-?-)eV 
\100kpc/\0.1 c)\\nG) 

on scales of R ~ 100 kpc, appropriate for the 
giant radio lobes in Cen A. Stochastic UHE 
proton acceleration in its lobes (as suggested 
in Hardcastle et al. 2009) would thus require 
Alfven speeds v A > 0.3 c. This seems diffi- 
cult to achieve (cf. also O'Sullvian et al. 2009). 
Thermal X-ray emission from the lobes Isobe 
et al. 2001 suggest (thermal) plasma densi- 
ties of n t h - (10~ 5 - 10~ 4 ) cirT 3 , implying va ~ 
0.003 c. Such values for « t h are consistent with 
independent estimates based on Faraday rota- 
tion measurements (Feain et al. 2009). Efficient 
UHECR acceleration in the lobes seems thus 



rather unlikely. 

(iii) Shear acceleration: If the flow profile is 
non-uniform across the jet, as in the case of a 
shear flow with u = u z (r)e z , then energetic par- 
ticles scattered across it, may be able to sample 
the flow difference du and thereby get acceler- 
ated (Jokipii & Morfill 1990; Rieger & Duffy 
2006). Like stochastic 2nd order Fermi, the av- 
erage energy gain is oc (du/c) 2 , although the 
physical origin is different (i.e., due to the sys- 
tematic, instead of the random motion of the 
scatterers). The velocity difference, that a par- 
ticles experiences, is du ~ (du z /dr)A, where A 
is the scattering mean free path. Again, this en- 
ergy change is acquired over t s ~ A/c, so that 



acc Ae/T s (du z (r)/dr) 2 A' w 

In contrast to eq. (4) and eq. (6), now f acc oc 
1/A. Thus, f acc becomes smaller as a particle 
increases its energy (so that its A oc p a , a > 0, 
becomes larger). Shear acceleration therefore 
preferentially picks up high-energy seed par- 
ticles for further energization, and acts more 
easily on particles of higher rigidity. Shocks, 
operating in the jet, could well provide the 
required seed particles (Rieger & Aharonian 
2009). Achievable particle energies are then 
constrained by the confinement condition (i.e., 
r gyro < width of the shear layer). The large- 
scale jet in Cen A has a projected length of 
about 4.5 kpc, and towards its end a width of 
about 1 kpc (Burns et al. 1983; Kraft et al. 
2002). For a characteristic B ~ lO" 4 /?; G on 
kpc-scale (cf. magnetic flux conservation) and 
a width of the shear comparable to the width of 
the jet, maximum energies 

E ~ ZeB(Ar) ~ 10 2 %Z eV (8) 

are possible. Shear acceleration might thus be 
able to boost energetic seed protons (e.g., pro- 
duced by shock acceleration) to energies be- 
yond E c . If the magnetic field gets amplified by 
internal shear (e.g., Urpin 2006), even bj ~ 1 
may be possible. Note that a shear dynamo 
could possibly explain why the magnetic field 
direction in Cen A seems to be almost parallel 
along its kpc jet (Hardcastle et al. 2003). 
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4. Conclusions 

The observed HE (Fermi) and VHE (H.E.S.S.) 
characteristics of Cen A ("core") suggest that 
the HE and VHE emission originate from dif- 
ferent regions. While the nuclear SED below a 
few GeV can be satisfactorily described with 
a conventional one-zone SSC model, an addi- 
tional contribution is required to account for 
the VHE emission. Non-thermal processes in 
the black hole-jet magnetosphere could offer a 
plausible explanation for the latter. 
Whether Cen A is indeed an UHECR source 
has observationally not been settled yet. From 
a theoretical point of view, efficient accelera- 
tion of protons to UHECR energies in Cen A 
remains challenging in the framework of most 
standard mechanisms. Observational evidence 
for UHE protons may therefore support the op- 
eration of an additional acceleration mecha- 
nism ("two-step") such as shear. The situation 
is much more relaxed for heavier elements like 
iron nuclei, which could most likely be directly 
accelerated (either by shocks or within the BH 
magnetosphere) to UHECR energies. Note that 
simultaneous operation of several mechanisms 
(with maximum energy not always linearly de- 
pendent on charge) also seems to constrain 
the potential to infer the UHECR composi- 
tion from the observed anisotropy (Lemoine & 
Waxman 2009). 
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DISCUSSION: 

VALENTI BOSCH-RAMON: Assuming 
ideal MHD in the jet launching region, how 
high can the particle Lorentz factor be? 

FRANK RIEGER: In ideal (single-fluid) 
MHD models, the azimuthal component of the 
magnetic field increases on approaching the 
light cylinder, and this makes centrifugal (test) 
particle acceleration inefficient. The details are 
dependent on the assumed rotation law and 
current distribution along the field line. It is, 
however, not always clear whether ideal MHD 
can be applied close to the BH (cf. Levinson & 
Rieger 2011). 



